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, Abstract. - Considering recent results revealing the existence of multi-scale rigid clusters of 

04 ' grains embedded in granular surface flows, i.e. flows down an erodible bed, we describe here 

the surface flows rheology through a non-local constitutive law. The predictions of the resulting 
model are compared quantitatively to experimental results: The model succeeds to account for 
the counter-intuitive shape of the velocity proflle observed in experiments, i. e. a velocity profile 
decreasing exponentially with depth in the static phase and remaining linear in the flowing layer 

ewith a velocity gradient independent of both the flowing layer thickness, the angle between the 
flow and the horizontal, and the coefficient of restitution of the grains. Moreover, the scalings 
, observed in rotating drums are recovered, at least for small rotating speed. 
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Introduction. - Granular media share properties with both usual liquids and solids: They 
' can sustain an inclined free surface without flowing but, when the angle exceeds a critical value 
identified since Coulomb § with some effective macroscopic friction angle <^m, an avalanche 
occurs. The motion has the peculiarity of being a surface flow. Two phases can be observed: A 
"solid" phase experiencing a creep motion where the averaged streamwise velocity decreases 
exponentially with the depth [^,^ , and a flowing phase exhibiting a linear velocity profile with 
a velocity gradient independent of the flowing layer thickness, the angle between the mean flow 
and the horizontal, and the coefficient of restitution of the beads [^|3| . Such profiles cannot be 
described using any conventional local and univocal constitutive laws of Continuum Mechanics 
for two main reasons: (i) The velocity gradient is found to be constant in the cascading layer 
whereas momentum balance implies that the shear stress increases linearly with depth and 
(ii) the velocity profile in the fiowing layer down an erodible bed differs significantly from the 
Bagnold-like velocity profile observed in dense flows down rough inclines ||]. Several models 
have been recently proposed to describe dense granular flows but, to our knowledge, none 
of them succeeds to account for the velocity profile measured experimentally in surface fiows. 
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A very recent experiment 0| has provided evidence of multi-scale rigid clusters embedded 
in the flow that may be responsible for non-local interactions in surface flows. Consequently, 
following the work of Mills et al. ^ originally devoted to dense flow down a rough incline, we 
propose here to describe the surface flow rheology through a non-local constitutive law. The 
resulting model is presented and discussed in the first section. In the second section, the model 
is applied to the description of ideal 2D flows at the surface of a packing of infinite depth. 
Among others, it accounts for the shape of the velocity profile observed experimentally. Finally, 
in the last section, we focus more specifically on steady surface flows in a slowly rotating 
large drum of finite thickness. The resulting model describes the experimental observations 
successfully. 

General formulation. - Let us consider surface flows of monodisperse spheres of diameter 
d and of density po in the gravity field g. Assuming that granular media deal with contin- 
uum models, the velocity and stress field should verify the momentum equation: pQi/ldtV + 
{v.'V)v} — pqvq + V.Z where v is the velocity of the fiuid particle, v the volume fraction 
and Z the stress tensor. In the following, variables are non-dimensionalized by length d, time 
\fdfg, and stress p^gd. Moreover, only steady, homogeneous 2D surface flows are considered: 
In the frame (cx, e^) tangent to the free surface (Fig. ^a), v(x^ z, t) = v(z)ex^ v(x^ z, t) — v(z) 
and the projection of the conservation equations gives: 

(6) £±fi = _^sin(? (1) 



dz dz 

A constitutive law should now be provided to relate these stress profiles to " observables" 
profiles such as velocity and volume fraction profiles. Correlated rigid clusters of grains em- 
bedded in surface flows have been experimentally evidenced . Their size has the peculiarity 
to be power-law distributed from the grain diameter to the flowing layer thickness. Conse- 
quently no characteristic correlation length scale - and thus no local stress/strain constitutive 
relationship - can be deflned. A non-local rheological law is thus proposed below. 

Following Mills et al. , we express the stress components as the sum of a motionless 
part sj^'* and an hydrodynamic part The motionless components are assumed to keep 

the relationship they got just before static equilibrium failure and the hydrodynamic parts 
are evaluated by considering the flow as transient solid chains embedded in a viscous fluid. 
This leads to (see Q for justifications and details): 

(a) Szj; = -|- siz' (5) "^xz — siz^ -|- S^z^ 

(c) ^=v^co^B (d) Si^i = -tan$Ei|^ (2) 

(e) =ao'-^J°a-iy)dy (/) l^i'J ^ + ^ iy)dy 

where ao is a constant coefficient related to the mean transient chains orientations 1^, C the 
correlation length of the chain network in a shear plane z = C*'', $ the internal friction angle 
of the material, ct" = rjdvx/dz the shear stress in the viscous fluid and vq the volume fraction 
of the immobile stack before failure that ranges from the random loose packed limit ~ 0.55 
to the random close packed limit v* ~ 0.64. Since the influence of the initial volume fraction 
Vq on the rheology is not investigated in this paper, we set = v* in the following. 

The internal friction angle $ is the sum of the particle/particle friction angle and of 
an angle (3 reflecting the entanglement of the grains |^ . According to the Reynolds dilantancy 
principle P vanishes when the volume fraction is smaller than the random loose packed 
value V* and increases with u when v > v* to be a maximum when v — v*. Consequently, 
$(i/) = $^ -I- (<i)J^^ - - I/*)/ (iy* - i^*)) where the maximum value $m of the internal 
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Fig. 1 - (a): Notation and sketch for the 2D surface flows The quantities R and 9 refer respectively 
to the flowing layer thickness and the mean angle of the flow, the transient solid objects are assumed 
to be separated by a mean distance constant in the whole packing. The frame {ex,ez) is chosen 
tangent to the free surface, (b) Rotating drum geometry. The radius of the drum and the solid rotation 
velocity are called respectively H and fi. The origin of {ex,ez) is now chosen to coincide with the 
center of the drum. For x — 0, dxV and dxVx can be assumed to vanish, but dxVz do not. The points 
Pi P2 denote respectively the upper and the lower boundary of the flowing layer. A point M of the 
free surface is defined through its curvilinear coordinate s whose origin is set at Pi. The surface S is 
normal to the free surface at M and bounded by the free surface and the drum boundary. 



friction angle can be identified with the angle of repose and f{x) is a dimensionless function 
increasing with repsect to x so that /(O) = and /(I) = 1. For sake of simplicity, / is identified 
with an Heaviside function to allow analytical solutions, i.e. : 

for V <Vi $ = $m, for V > Vi $ = $m (3) 

Expressions for the viscosity are usually obtained from the Enskog equation It reads 
77 oc g(y)\fT where g{v) is the pair correlation function expected to diverge as (y*^ — v)^^ in 
the high density limit [|l2[ and -v/T is the granular temperature identified with the rms velocity 
fluctuations. For dilute flows, the temperature profile can be determined by a local balance 
between the viscous heating and dissipation through inelastic binary collisions But in 
dense flows, the kinetic energy of a collision is dissipated within the entire packing through 
multiple collisions and the concept of a local energy balance becomes irrelevant P, |^, . 
Moreover, the mean free path is very small and the experimental noise is expected to be 
sufficient to drive collisions jlj]. For all these reasons, we choose here to define an effective 
temperature Te// constant in the whole packing, resulting from the global energy balance 
within the whole system. The viscosity can then be written as: 

= '7o5('^)-T^ where g{iy) = — (4) 



Finally, from Eqs. |l] and H, one can relate v to a"" via: 



aotan<i> „ 
Ccos0 



(5) 



Canonical case: Flows down to the surface of an infinite packing. - Let us consider now 
the granular flow at the surface of a packing of infinite depth and infinite width. From Eqs. ^ 
and H, the viscous stress should obey: 
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where the functions A and F are given respectively by: F{9i,92) — cos 6*1 (tan 6*1 — tan6'2) 
and A(0i,02) = C/(l + ao tan^?i tan6'2)- For sake of simpHcity, 9 variations are assumed to 
be smah around and the different trigonometric functions are expanded in the relevant 
order of 6* — ^m- The expressions of F and A become: F{9i,92) — {9i — 9-2)1 cos$m and 
A(0i, 02) ^ C/(l + «o tan2 $m) ^ C*^ 

For a surface flow, <J^{z) is expected to vanish for infinite depth. From Eq. ^, Eq. ^ and 
Eq. one deduces sucessively that i/ ^ i^*, $ ^ $m and F{9,^m) when z —oo. 
Consequently, permanent homogeneous surface flows are available in infinite packing only for 
9 = <i>M- When z is gradually increased from — oo, first $ = <I>m (Eq. ||), ct" increases (Eq. ^ 
and 1/ decreases (Eq. ||). Then, at a given depth identified with the cascading layer /erodible 
bed interface, v — h'l , which switches the value of $ = $m to $ = (Eq. ||). Consequently, 
Eq. 1^ can be written: 

for ze[~^,^R] ^-la-^O 

for z e [-R, 0] ia" - -i^:F{^m, $™) ^ ' 

where R is the flowing layer thickness. Let us note that, in this approach, the presence of 
the two phases, the "static" phase and the flowing layer, comes from the variation of $ with 
respect to ly. Then, as , S^^^ and consequently a" (Eq. |^) should vanish at the free surface, 

and as S^^' , and consequently cr" (Eq. H) should vary continuously at the interface , one 
gets: 

for z e [-(X3, -R] (T^'(z) = ;.:AF($m, $™)(e^/^ - l)e'/^ . . 

for z e hi?, 0] a"{z) = iy*XF{<S>M, $m)(l - e^A) ^ ' 

In the following, one assumes A <C i? and does not consider what happens in the vicinity of 
the free surface — A < z < 0. The tr'^ profile is then given by: 



for z e [-00, -R] a^iz) ~ j/*AF($m, *m)e^/^e^/^ 
for z e hi?, -A] a-{z) ~ z/*AF($m, $m) 

From Eqs. |^ and ^ the volume fraction profile can be deduced: 



for ze[-oo,-R] Kz) = ^*(l-(l-A)£^-i:^e«/V/^ 



C ' sin $4 

\\ F{<S>M,'l 
Q ' sin 



for z e h^, -A] i^iz) = ( 1 - (1 - 



(9) 



(10) 



The volume fraction is thus predicted to be different, but constant in the two phases except 
in a narrow layer at the interface where the two values are connected through a matching 
exponential function of characteristic length A. From Eqs. and |lo| the velocity profile can 
be deduced: 



for z € [-00, -R] v^{z) = A^oe2«/^e2^/^ 
for 2 6 hi?, -A] v4z) ^ io{z + R + ^) 

where the velocity gradient 70 is given by: 



(11) 



^^F(£m^ A (12) 
r,o sin$M ^ C ^ ^ 
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which captures the main experimental observations, namely a velocity profile decreasing ex- 
ponentially with depth in the static phase |^,^ and remaining linear in the flowing layer with 
a velocity gradient independent of i? [^,^- 

Let us note finally that the use of a generic viscous law a"" — rjo^dvx/dz)" / modifies 
the shape of the velocity profile only through 70 and the characteristic length of the exponential 
term that become respectively 70 ^ {iyf^'^F{^M, $m)''+^A(l - A/C)'^/(77o sin^ $a/)}^/" and 
A — > a A/ {P + 1). This indicates that the precise expression of the viscosity, and in particular 
the assumption of a constant temperature, has no influence on the shape of the velocity profile. 

Steady surface flows in a large rotating drum. - Let us now investigate permanent surface 
flows at the center of an half-flUed rotating drum of radius H (cf . Fig. |l|b) . The origin of the 
frame {ex,ez) is now chosen at the center of the drum with (resp. e^) is tangent (resp. 
normal) to the free surface (see fig. |l]b). The rotating speed Q is supposed to be small so that 
only first order terms in are considered. The inertial stresses scale as ft^ and can thus be 
neglected in Eq. |l|. In such geometry, the flow cannot be considered as homogeneous anymore. 
d^v and d^Vx vanish for x — 0, but d^v^ do not: At the drum boundary, v is given by the 
solid rotation of the drum which leads to dxVz{x = 0;z = —H) = —VI. More generally, one 
can write Vz{5x]z) = —flSxf{z) where / is a dimensionless decreasing function of z so that 
f{—H) — 1 and /(O) = 0. This introduces new terms = = —rjyV.v = —rj^dzVz in 
the stress tensor representing the dissipation generated by dilation where rji, is the volume 
viscosity and an unknown function of v. The conservation equation previously given by Eq. ^ 
becomes dx^^x + dS^z/dz — —I'smO and dEzz/dz — vcosO. Then, cr" should obey: 

^ - = -i^lFnie, V, z) (13) 

where Fa{9,^,Vt,v,z) = F{0,^) + r]^{iy)flf'{z)/v*. As (j''{-H) vanishes, we deduce that 
v{-H) = V* (Eq. |), $(-iJ) = $M (Eq. |) and, assuming that iJ > i?, Fn{e,<^M, Vt, i/*, -H) ~ 
(Eq.H), I.e. : 

= $M + pn where p = cos $m | /'(-i?) | (14) 

which reproduce the scaling relations 0{fl) observed experimentally when Q is small enough 
to neglect inertial effects 

The forms of f{z) and are now needed to evaluate the velocity profile. For sake 

of simplicity, one chooses here the simplest form compatible with the boundary conditions 
/(O) = and f{—H) — 1, i.e. f{z) = —z/H. For the same reason, rjs is also set constant in 
the static phase. Then, tr" obeys: 



for ze[-H,-R] ^-{a- = -iy:Fie,'i>M)+vM/H = 

for z e [-R, 0] {a- = ~'^*F{0, $,„) + ^M/H = -i^*F{^m,<^>„,) ^ ' 

And the volume fraction profile (resp. the volume fraction profile) are given by Eqs. ^ 
and ^ (resp. Eq. p^ ). The viscous relationship now reads ct" = riog{h'){dvx/dz + dxVz) with 
dxVz — ilz/H and, assuming that » i? » A, the shape of the velocity profile is finally 
given by: 

for ze[-H. -R] vx (z) ~ |7oe2«/^e2z/A _ ^ (^2 ^ ^2) 

for ^e[-i?,-A] Vx{z):^'yoiz + R+^)-^{z' + H^) ^ ' 
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Fig. 2 - Comparison between theory and experiment. (a): Variation of the angle 6 (in radian) of the 
mean flow versus the rotation speed Q (non-dimensionalized by ^ gjd). Circles correspond to the 
experimental results of Q obtained in a rotating drum of diameter H — 45 cm and of thickness 7 mm, 
half-filled with steel beads of diameter d = 3 mm. The critical angle $Af has been measured to be 
$A/ = 27.5°. The full line is the theoretical result given by Eq. ^ with a fit parameter p = 32.6. (b): 
Velocity profile v^^iz) in non-dimensionalized Units. Full lines corresponds to theoretical profiles given 
by Eq. Experimental data are fitted by changing the parameter 70 for the entire set of profiles 
and by evaluating R from Eq. |l^ for each profile. The fitted correlation length and velocity gradient 
are found to be A = 3.62 and 70 = 0.63 respectively. 



And consequently, for <C 70 as in experiments [y|3|, the solid rotation of the packing 
does not modify significantly the shape of the velocity profile. The relationship between the 
flowing layer thickness R and can now be deduced: For a steady flow, at each point M 
of the free surface, the flow rate (5(M) through the surface S'(M) normal to the free surface 
at M and bounded by the free surface and the drum boundary (Fig. |^b) should obey the 
mass conservation equation: dsQ — where s is the curvilinear coordinate of the point M. 
As Q vanishes at the upper and lower boundaries of the flowing layer (points Pi and P2 
of Fig. |l|b), Q vanishes everywhere, in particular at the center of the drum. One thus gets: 
Q = v,j:{z)dz = 0, and, assuming X ^ R <ti H: 



1/2 



which is in agreement with the experimental scaling reported in [ i , |3|, [l^ . 

In Fig. ^, the model predictions are compared to the data of |3| measured experimentally 
in a rotating drum of diameter H = A5 cm and of thickness 7 mm, half-flUed with steel beads 
of diameter d — 3 mm with a particle/particle friction angle = 10°. Moreover, the angle of 
repose has been measured in the experiment: <I>m = 27.5°. Figure ^ja presents the variation of 
with measured experimentally in |^ and calculated using the theoretical scaling (Eq. |lj) . 
The flt parameter is found to be p = 32.6. Figure ||b presents the comparison between the 
velocity profile given by Eq. |l^ and the ones measured experimentally in Q . The last ones are 
fitted by setting H = 75, by evaluating R using Eq. [T^ for each profile, and by changing 7-0 and 
A for the entire set of profiles. The agreement is remarkable. The value of the fit parameters 
are found to be 70 = 0.63 and A = 3.62. 

Conclusion. - In this letter, we have proposed a non-local diphasic model somewhat 
derived from M to describe surface flows in granular packing. This model accounts for the 
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shape of the velocity profile experimentally observed, which no theories allowed for until now. 
Then, we have applied the model to the description of steady surface flows in a rotating drum. 
All the experimentally observed scaling were recovered, at least for rotating speed slow enough 
to let inertial effect be negligible. Our theoretical predictions agree quantitatively with the 
experimental results of 

A main limitation of the model comes from the assumption that the two lengths charac- 
terizing the granular mesostructure, i.e. the length ^ of the transient chains (infinite here) and 
the average distance A between the chains in a shear plane, are fixed throughout the entire 
packing. Assuming connectivity transition in a finite size system, one can relate these two 
lengths to the volume fraction and the size of the media. This allows to extend the present 
non-local description to other geometries such as flows down rough inclines and shear 
Couette flows . This will be explored in a forthcoming publication . 

* * * 

We acknowledge discussions with E. Bertin, F. Chevoir, F. Daviaud, O. Dauchot and L. 
Laurent. We also thank Dan Qvale for the careful reading of this manuscript. 
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